RADIAL CORRELATIONS BETWEEN TWO QUARKS 
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In nuclear many-body problems the short-range correlation between two nucleons is 
well described by the corresponding correlation in the iu;o-body problem. There- 
fore, as a first step in any attempt at an analogous description of many-quark 
systems, it is necessary to know the two-quark correlation. With this in mind, 
we study the light quark distribution in a heavy-light meson with a static heavy 
quark. The charge and matter radial distributions of these heavy-light mesons are 
measured on a lattice with a light quark mass about that of the strange quark. 
Both distributions can be well fitted upto r 0.7 fm with the exponential form 
iii?(r), where Wi{r) = Aexp(— r/r^). For the charge(c) and matter(m) distributions 
rc ~ 0.32(2)fm and rm ~ 0.24(2)fm. We also discuss the normalisation of the total 
charge (defined to be unity in the continuum limit) and matter integrated over all 
space, finding 1.30(5) and 0.4(1) respectively for a lattice spacing 0.17 fm. 

1 Correlations between two nucleons 

In nuclear many-body problems it is well known that the short-range correla- 
tions between pairs of nucleons are often well described by the corresponding 
correlations of the two-body problem. For example: 

• Few nucleon systems. The NN correlations for the deuteron, '^He 
and *He, when normalised to each other at their maxima, are essentially 
indistinguishable upto about 2 fm. Only for larger internucleon distances 
do the binding energy differences play a role El. 

• Many-nucleon systems. Often variational methods are used and a 
typical trial wavefunction is taken to be of the form 

|*T [I^^<JF,,] |*(SheU Model) >, (1) 

where ^ (Shell Model) is simply a product (antisymmetrised) of single 
nucleon wavefunctions - encoding the long range correlations. However, 
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the short range physics between the NN pairs is in the 



F — 



(2) 



where = [1, ai.aj, ...] {^[1, Ti.Tj] are the various spin and isospin op- 
erators for two nucleons. The variation then yields radial correlatio 
Umi^ij) that are ~ 2-body problem correlations - see for example Ref. 

• Nuclear Matter An extreme use in many-body systems fii the two- 
body correlation is the Moszkowski-Scott separation methodd in nuclear 
matter. There the many-body wavefunction is taken to be of the form 

*(Many - Body) = ^(Two - Body) r < 

= sin(fcr) r > rc, 

where Tc is the value of r at which ^'(Two — Body) and sin(fcr) and their 
derivatives equal each other. 

The conclusion is that, in many-nucleon systems, the tvifo-nucleon 
correlation from the two-body problem plays an important role. 

The question can then be asked whether similar effects arise in miilti- quark 
systems. This may help in the understanding of such systems. At present and 
in the forseeable future, since lattice QCD calculations are restricted to very 
few quarks (upto four), some model that encodes the results of these few- quark 
calculations is needed in order to proceed to systems containing more quarks. 

Earlier attempts by the present authors have concentrated on understand- 
ing the energies of four-quark systenia in terms of interquark potentials. These 
involved four heavy quarks {Q'^Q^p and the interaction of two heavy- light 
mesons (Q^q^) H. In neither case was it possible to understand the four-quark 
energies in terms of simply two-quark potentials - with the result indicating 
that a four-quark form-factor was necessary. This conclusion has also been 
supported by colour field measurements Here, as a first step, the emphasis 
is on extracting the radial correlation between the two quarks of a heavy-light 
meson. Hopefully, this additional information on the Qq system will enable 
a less ambiguous model of this system to be made. At present, if through 
an operator 0{r), a transition between two two-quark states is calculated as 
{tpi{r)\0{r)\ip2{r)) , then the ipi are first calculated as a byproduct of some dif- 
ferential equation containing some interquark potential - say V{r) ~ a/r-\-br. 
However, both the form of V{r) and the equation into which it is inserted 
are not unique. The latter could range from a non-relativistic Schroedinger 
equation, if the quarks are heavy enough, through a series of equations that 
incorporate relativistic effects to varying degrees. 
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2 Heavy-light mesons [Qq) 

A study of heavy-light mesons is not just for the academic interest dis- 
cussed above, since they are realised in nature. The best examples are the 
5(5.28 GeV) and i?s(5.37 GeV) mesons, which have quark structures (bu) 
and(6s). Since the 6, s and u-quarks have masses of about 4.2, 0.1 and 
0.001 GeV, we see that the B and Bg are indeed heavy-light mesons and 
can be thought of as the "Hydrogen atom" for quark systems. Currently the 
i?-mesons are of particular interest, since they are expected to lead to a bet- 
ter understanding oLCP violation and for this reason are being generated at 
so-called B-factoriesQ. 

More details of the following sections can be found in Ref.l3. 



2.1 Energies of Heavy-light mesons 

The basic quantity for evaluating the energies of heavy-light mesons is the 
2-point correlation function - see Fig. |l| a). 

C(2, T) = ({/^(x, i, T)P(x, t + T, t)}, (4) 

where U'^ {x,t,T) is the heavy(infinite mass)-quark propagator and P{x,t + 
T, t) the light-quark propagator. The (...) means that C(2) has been averaged 
over the whole lattice. Since the C(2) decay as exp{—EoT), where Eq is the 
energy of the ground state, we get 

These energies were calculated in Ref.i for the states 
Lj = Si/2, Pi/2, P3/2, D3/2, -D5/2, F 

2.2 Heavy-Light radial correlations 

Here the basic quantity is 3-point correlation function ~ see Fig. ^ b). 

C(3,-t2, h, r) = (C/'3(x,-t2, ti)Pi(x,ii;r,0)e(r)P2(r,0;x,-i2)), (6) 

where the Pi. 2 are the light anti-quark propagators that go from the Q at time 
ti to the point r at i = and then return to Q at time — i2- The probe 8(r) 
at r is here considered to have two forms i) = 74 for measuring the charge 
distribution of the q and ii) = 1 for measuring the matter distribution. 
Knowing the C(3) then the radial distributions are given by 
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2.3 Lattice Parameters and Refinements 

The essential lattice parameters are: 

1) Lattice size 16^ x 24 

2) Quark-gluon coupling /3 = 5.7 giviag a lattice spacing « 0.17fm 

3) Light quark mass « mg - see Ref. llJ as evidence for this. 

4) Quenched Approximation i.e. quark-antiquark pairs are not included. 
Eqs. 1^ and |^ for the ground state energy(i?o) and correlations F[C(0, T, R)] 
are based on a single state constructed from the basic lattice. However, 
this can be extended by modifying this original lattice to give other states 
(i,j,k...). In this way the 2- and 3-point correlation functions can become 
matrices Cij{2), Cij{3, r). Their diagonalisation - by reducing the contami- 
nation of higher states on the ground state - then gives improved estimates of 
i?o and F[C{Q, T, R)]. The mechanism for generating these addition states is 
fuzzing. This replaces a link on the lattice by a combination of neighbouring 
links to give a fuzzed link - with projection to SUB implied ~ as 

[A fuzzed link] = /p- [Straight link]-)- [Sum of 4 spatial U-bends]. 

This cycle can be repeat many times. Here, in addition to the basic 
state(L), two new states Fl and F2 are constructed using 2 and 6 cycles with 
fp = 2.5. respectively. 



3 Results 



The lattice data {dj (2, T)) and (dj (3, T, r)) can be analysed in several ways. 
Here two are presented: i) Visual and ii) Fitting with exponentials. 



3. 1 Visual 

Fig. I shows from Eq. ^ the ratio (C(3, i?)}/(C(2)) for values of r upto 5 
lattice spacings i.e. about O.Sfm. As T ^ oo, these clearly show plateaux, 
which give directly the desired charge or matter density. 



3.2 Fitting with exponentials 

A more precise method is to first fit the (Ci.j(2, T)) by the approximate ex- 
pression 

3 

a, (2, T)^Y. < exp(-i?„r)«;, (8) 

where the i, j refer to the states L,F1,F2. This results in the energy eigenvsilues 
Ea and their eigenvectors . These are the energies quoted in Ref. It 
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Table 1. The parametrization of the charge and matter densities as ui^(r), where Wm(r) = 
Am exp[-T-/r„i]. 





A, 


Ti (fm) 


RMS Radius (fm) 


Charge Density 
Matter Density 


0.26(1) 
0.29(1) 


0.32(2) 
0.24(2) 


0.55(3) 
0.42(3) 



is found that only data with T > 4 gives /dof k, 1. This means that 54 
data points arc fitted with 12 parameters. 

Given the Ea and vf the Cij(3, T, R) are then fitted by the approximate 
expression 

3 3 

C',,(3,r,r) = ^ ^<exp[-£;„ti]a;"'5(r)exp[-i;;3i2]<, (9) 

a=l 13=1 

where the x°'^{r) are varied giving directly the desired charge or matter den- 
sity. This needs ii + i2 > 8 to get /dof k, 1 when fitting 18 data points 
with 3 parameters. 

3.3 Analytic forms for the charge and matter densities 

For the above charge {i=c) and matter {i — m) densities, it is of interest to 
express Xi{r) = wf{r) in the form Wi{r) = yl^ exp(— r/r^), where the parame- 



ters ri acd Ai are given for a global fit in Table 3.3. There, as in Fig. it is 
clearly seen that the charge density has a significantly longer range than that 
of the matter density. 

The fourier transforms of these densities lead to vector(z=ti = c) and 
scalar(i = s = to) form factors that are usually expressed as [F^ sio^)] 
{q^ + Ml_g)^^ . These are appropriate for the long-range part of w1{r) and 
lead to My—Q.%{l) and Mg = 1.3(l)GeV and should be compared—with the 
the direct calculation of the ss vector and scalar mesons in Refs. lid and O. 
The latter only involve 2-point correlation functions and so are more precise 
leading to M„=0.944(2) and Ms = 1.61(6)GeV. 

3.4 Sum-rules 

The above measures the densities at a few definite values of r. However, it is 
of interest to consider the sum-rules that sum over all values of r i.e. 

F^-^[C{2,)] = ^^Ic^^^^^^^ - (10) 
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For the charge and matter densities this yields Jc=1.30(5) and /,„=0.4(1) 
respectively. On the other hand, integrating the above analytic expressions 
for Wi(r) gives 1-5(1) and 0.7(1). These estimates are less reliable than the 
direct lattice sum, since much of the contributions are from small values of r 
where lattice artifacts enter. 

By charge conservation we should get /c=l in the continuum limit - 
the quark charge having been defined as unity. However, it is known that 
renormalisation effects enter due to the finite lattice spacing. In Ref. t3, 
for lattice parameters different to those used here, an overall renormalisation 
factor {F^) of « Q_S is found for the vector (charge) vertex. Similar reductions 
are found in Ref. i3 for the axial vector operator. The inverse of our value of 
Jc=1.30(5) could then be interpreted as a non-perturbative estimate of as 
0.77(3). In any case the conclusion is that our /c=1.30(5) is consistent with 
charge conservation for the continuum. 

3.5 Dirac amplitude interpretation 

In a relativistic quantum mechanical picture the above charge and matter 
distributions can be interpreted as + |/P and |gp — where g, f are the 
upper and lower components in the solution of the Dirac equation. Possibly 
this will now serve as a way to remove some of the ambiguities refered to at 
the end of Section |l] by ensuring that the solutions of the Dirac equation agree 
not only with the lattice or empirical energies but also the forms of !<?(?') P 
and |/(r)|2. 

Conclusions and Future 

• The S'i/2-wave charge and matter densities can be measured quite reliably 
out to « G.7fm. 

• Off-axis measurements can improve these results. In particular, it is in- 
teresting to see whether or not the charge density for r = 5 is indeed lower 
than the simple exponential drop off in Fig. ^ This would indicate that 
the confining potential is playing a role. Hopefully, the better statistics 
at the off-axis point (3, 4) would clarify this point. 

• The Pi/2; -P3/2J -D3/2, -D5/2, ... densities are now being measured. 

• Understand the densities phenomenologically using the Dirac equation. 

• Measure correlations in the {Q^(f) system and check their form against 
the above Qg-correlations. 
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• Use larger /3 and lattices to get nearer the continuum limit. 

• Replace the present quenched lattices by unquenched ones. 

• Consider other operators, e.g. Pseudo-vector (7^75) for B*Btt couphng. 
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Figure 1. a) A two-point correlation function, b) A three-point correlation function. 
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Figure 2. The ratio (C(3, r))/(C(2)> for 
r = 0,...,5. 
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